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ABSTRACT 

The high velocity dispersion compact cloud CO—0.30—0.07 is a peculiar molecular 
clump discovered in the central moleculr zone of the Milky Way, which is characterized 
by its extremely broad velocity emissions 145 kms -1 ) despite the absence of internal 
energy sources. We present new interferometric maps of the cloud in multiple molecular 
lines in frequency ranges of 265-269 GHz and 276-280 GHz obtained using the Sumb- 
millimeter Array, along with the single-dish images previously obtained with the ASTE 
10-m telescope. The data show that the characteristic broad velocity emissions are 
predominantly confined in two parallel ridges running through the cloud center. The 
central ridges are tightly anti-correlated with each other in both space and velocity, 
thereby sharply dividing the entire cloud into two distinct velocity components (+15 
kms -1 and +55 kms -1 ). This morphology is consistent with a model in which the 
two velocity components collide with a relative velocity of 40 kms -1 at the interface 
defined by the central ridges, although an alternative explanation with a highly inclined 
expanding-ring model is yet to be fully invalidated. We have also unexpectedly detected 
several compact clumps (< 0.1 pc in radius) likely formed by shock compression. The 
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clumps have several features in common with typical star-forming clouds: high densi¬ 
ties (10 6 - 5-7 - 5 cm -3 ), rich abundances of hot-core-type molecular species, and relatively 
narrow velocity widths apparently decoupled from the furious turbulence dominating 
the cloud. The cloud CO—0.30—0.07 is possibly at an early phase of star formation 
activity triggered by the shock impact. 

Subject headings: Galaxy: center, ISM: clouds, stars: formation 


1. INTRODUCTION 


Widespread strong turbulence is an outstanding characteristic of the molecular clouds (MCs) 
in the central molecular zone (CMZ) of the Milky Way. Their velocity dispersions lie about a factor 
of 5 a bove the size-velocity width relationship for the Galactic dis k, irrespective of choice of tracer 
lines f Mivazaki fc Tsuboi 2000l : Oka et al. 2001 : Shettv et al. 20121 ). The virial parameters (<a v ir) of 
the CMZ clouds are systematically ^10 times those in t he Galactic disk, indicating that they are no t 
self-gravitating but likely bound by external pressure (IMivazaki fe Tsuboill2000l : lOka et al.ll200ll ). 
It is suggested by many authors that the high degree of turbulence in the CMZ may affect the star 
formation process there. Primarily, the strong turbulence suppresses the star formation rate (SFR) 
because the gravitational collapse of dense clumps in the MCs is slowed or inhibited by the high 
turbulent pressure. Several recent works suggest that the SFR within the CMZ may be suppressed 
by thi s effect to a value 1-2 orders of magnitude lo wer than the value expected from its dense gas 
mass (ILongmore et al.ll2012bl : iKruijssen et al.ll20141 h although other effects such as the tidal shear, 
strong m agnetic field, and enhanced cosmic-ray flux m ay also play significant roles in determining 
the SFR ( Yusef-Zadeh et al. 2007 : Kruijssen et al. 20141 ). Meanwhile, turbulence may also facilitate 
the formation of massive stars and stellar clusters in sites of cloud-cloud collision (IHabe fc Ohta 
1992l : llnoue fe Fukuill2013l ). This collision-induced s tar formation is a pop ular explanation for the 
highly active star formation in the Sgr B2 complex (jHasegawa et al.1119941 ). and there is increasing 
evidence suggesting that cloud-cloud collisions play a key ro le in cluster formation in the CMZ 


(Hig uchi et al.l 120141 ) as well as in the Galactic disk region (IFurukawa et al .1 120091 : llliguchi et al 

20Kj). 


The properties and origin of turbulence in the CMZ are, however, yet to be fully understood. 
Shettv et al.l (2012) show that the size-linewidth relationship holds both within and among the 
CMZ clouds over an order of magnitude in the spatial scale. This indicates that the turbulence 
driving in the CMZ is predominantly a large scale one, which is consistent with the idea that the 


turbulence is maintained by interactions with frequent supernova explosions (Shettv et al.ll2012l ) 
or with the bar potential in the inner Galaxy (jRx)driguez-Fernandez et al.ll2006l ). Meanw hile, high 


spatial resolution surveys in millimeter and submillimeter carbon monoxide (CO) lines (Oka et al 


1999, 20121 ) detected a considerable number of small clumps with extremely broad velocity widths 
inconsistent with the size-velocity width relationship; their velocity widths are 40-120 kms -1 
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in full width at zero-intensity (FWZI) in spite of their small spatial sizes of approximately 1-3 
pc in radius, indic ating e nhancement in velocity widths by a factor of 4 from the size-linewidth 
relationship (Tanaka et al. 2014). These clumps are distinguished by their spatial compactness 


Bania 


from the larger-scale_broad-line structures of a few 10 pc in sizes that were identified by 
ljl.977h and ll is/tl fepQfih at lar ger Galacti c longitudes (|Z| > 1 °). The recent CO J— 3-2 survey 
with the ASTE 10-m telescope ( Oka et al. 20121 ) identified 70 such high velocity dispersion compact 


clouds (HVCCs; Oka et al. 2007) and HVCC-like features, most of which were not associated with 
known energy sources for driving such small scale turbulence. 

Recently, we reported results of the HCN J=4-3 mappings obtained with the ASTE 10 -m 
telescope toward the Sgr C complex, which include one of the most prominent HVCCs in the 
CMZ, CO—0.30—0.07 (Tanaka et al . 2014). We showed that the cloud had a velocity width of 
120 kms -1 comprising a pair of broad velocity lobes with a velocity width of 60 kms -1 each, 
which were anti-correlated with each other in both space and velocity. The kinetic energy of the 
entire system was estimated to be 10 49 ergs, which is comparable to those of the clouds interacting 
with supernova remnants (SNRs). We proposed several hypotheses for the origin of the HVCC: 
cloud-cloud collision, hyper-energetic molecular outflow, a compact expanding shell, and a rotating 
ring/disk, although the spatial resolution of the telescope (24" = 1 pc) was insufficient to distinguish 
these possible cases. 

This paper reports the results of new observations toward CO—0.30—0.07 performed with the 
Submillimeter Array (SMA). We present interferometric images of the cloud in submillimeter lines 
of HCN, HCO + , N 2 H+, CH 3 OH, and several other hot-core-type molecules probing the dense molec¬ 
ular gas (n > 10 5 cm -3 ) in a wide range of physical and chemical conditions. The HCN, HCO + , 
and N 2 H + lines are commonly used high density tracers, which are sensitive to diff eren t environ¬ 
ments; N 2 H + tends to selectively trace quiescent, cold cores (e.g. Tatematsu et al. 2008), whereas 
HCN and HCO + abundances are less environmentally depe n dent and hence they are often useful 
in detecting shock-compressed gas (e.g. Tanaka et al. 2007, 2014). Methanol is one of the most 
abundant interstellar organic molecules with enhanced abundance in star-formin g cores , though 


other enhancement mechani sms such as cosmic-ray-in duced photo-desorption (lYusef-Zadeh et al 


2013) and ejection by shock (IReauena-Torres et al . 2006) are also suggested to be important in the 
CMZ. With these data, we investigate the structure of the shocked gas in the cloud in detail and 
thereby examine whether it shows systematic motion such as collimated outflow, expansion, and 
rotation, or its large velocity dispersion is caused by random turbulent motion. We also report the 
unexpected discovery of compact clumps toward the boundary between the two velocity compo¬ 
nents of the cloud, which may indicate an early phase of star formation activity under the highly 
turbulent envi ronment. In this paper, we adopt 8.3 kpc for the distance to the Galactic center 
(Gillessen et al. 2009). 
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2. OBSERVATIONS 
2.1. SMA Observations 


The observations were performed on May 12 th, 2013, by using the Submillimeter Array (SMA) 
in the sub-compact configuration providing baseline lengths ranging from 9.5 m to 25 m. We 
conducted a hexagonal 7-point Nyquist-sampled mosaic to cover the 52" radius region centered 
at (£* 2000 ,^ 2000 ) = (17 h 45 m 08 s .0, —29° 13'44".8) or (l,b) = (-0°.3053, -0°.0622) in the Galactic 
coordinates. The 230 GHz receivers were tuned to observe frequency ranges from 264.77 GHz to 
268.75 GHz (LSB) and from 276.77 GHz to 280.75 GHz (USB) simultaneously. The correlator 
configuration was chosen to provide 4 GHz frequency coverage in each sideband and 812.5 kHz 
spectral resolution. The HCN <7= 3-2, HCO + <7= 3-2, and CH 3 OH Jr— 6 i -52 E, 52 ~ 4 i E, and 9o- 
81 E lines were observed in the LSB, and the N 2 H + <7= 3-2, 13 CS <7=6-5, OCS <7=23-22, 0 -H 2 CS 
Jk&,k c — 8 i 5 7 - 7 i 5 6 , and CH 3 OH <7k=9_i-8oE lines were in the USB. We list the target lines in 
Table [U along with their rest frequencies and upper state energies. 


We observed Neptune, 3C279, and nrao530 for the flux, bandpass, and gain calibration mea¬ 
surements, respectively. We employed the CASA package developed for the NRAO for flagging bad 
data, applying calibrations, and creating images. We applied Briggs weighting for the imaging, 
which provided a 8". 2 x 5".2 synthesized beam for HCN <7= 3-2. The beam position angle is — 0°.5 
in the equatorial coordinates or — 59°.1 in the Galactic coordinates. Phase-only self-calibration was 
applied to the quasi-continuum data in each sideband, the calibration table of which was trans¬ 
ferred to the line data. The continuum image was constructed by integrating the emissions of the 
line-free channels in the total 8 GHz bandwidth. The total on-source integration time was 25 min 
per pointing, and the resultant r.m.s. noise level of the final maps is 0.22 Jybeam -1 per velocity 
channel for the USB, 0.27 Jybeam -1 for the LSB, and 3.9 mJybeam -1 for the continuum data. 
Correction for the primary beam was not applied. The maps were originally made in the equatorial 
coordinate system, and subsequently we converted them into the Galactic coordinate system for 
convenience of com parison with the single-dish HCN <7=4-3 data obtained with the ASTE 10-m 
telescope (Tanaka et al. 2014). 


2.2. ASTE Observations 


We also use t he published data of HCN <7=4-3 and 13 CO <7=3-2 obtained with the ASTE 10 -m 
telescope (Tanaka et al. 2014). The maps were obtained with 24" spatial resolution and 2 kms -1 
velocity resolution by employing the On-The-Fly (OTF) scan mode. The r.m.s. noise levels in the 
main-beam temperature scale are 0.13 K and 0.21 K for the HCN and 13 CO data, respectively. 
Further details of the ASTE data and observations are described in lTanaka et al.l (2014). 
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3. RESULTS 

3.1. Overall Structure: Central Clump and Extended Turbulent Component 

Figure [T] shows the peak flux densities of the observed lines and continuum flux density, where 
the peak flux densities were calculated from spectra binned over 10 velocity channels to improve 
S/N ratios. Two spatial components can be recognized in the maps. One is the bright compact 
clump toward the field center, which is commonly seen in all the observed lines and the continuum 
but more clearly in the maps of the optically thinner group (i.e., the continuum and the lines other 
than HCN and HCO + ). This central clump is point-source-like in the methanol, OCS, 13 CS, and 
0 -H 2 CS maps, whereas its N 2 EU and continuum images have sizes significantly larger than the 
beam size. The other component is the clumpy, spatially extended component within a radius of 
approximately 30" around the central clump. This extended component is dominant in the HCN 
and HCO + maps and is recognizable in the maps of N 2 H + and CH 3 OH 52 ~ 4 i E and 9_i-8o E as 
well as the continuum. The mass distribution is better represented by the optically thin group, 
indicating that the major fraction of the cloud mass is confined in the central clump. 

The contrast between the central clump and the extended component is more striking in Figure 
El which shows the posit ion-velocity (P- V) diagrams along the Galactic latitude axis for the SMA 
data of the HCN and N 2 H + lines averaged over the full latitude range. The extended component in 
the HCN map shows a highly turbulent velocity structure with emission spanning in a tlsr range 
of 145 kms -1 from —35 kms -1 to +110 kms -1 , which is wider than that previously measured with 
our single-dish observations. In contrast, the central clump in the N 2 H + image is strictly limited 
to a narrow velocity range from —10 kms -1 to +10 kms -1 , apparently decoupled from the outer 
highly turbulent flow. The HCN image has a peak toward the central clump at a velocity slightly 
blue-shifted by 5 kms -1 presumably because of self-absorption. We also show the single-dish HCN 
J=4-3 and 13 CO J— 3-2 data in Figure El in which only the extended component is visible; they 
do not have a narrow peak at the velocity corresponding to the central clump. 


Table 1. Molecular Lines Observed with the SMA 


Molecule 

Transition 

Frequency 

(GHz) 

E u /kB 

(K) 

sideband 

USB/LSB 

HCN 

3-2 

265.886 

25.5 

L 

HCO+ 

3-2 

267.558 

25.7 

L 

n 2 h+ 

3-2 

279.512 

26.8 

U 

13 cs 

6-5 

277.455 

23.1 

U 

0 -H 2 CS 

8l,7-7l,6 

278.888 

73.4 

u 

OCS 

23-22 

279.685 

161.1 

u 

CH 30 H 

6i-5 2 E 

265.290 

61.9 

L 


5 2 -4i E 

266.838 

49.2 

L 


9o-8i E 

267.403 

109.6 

L 


9-i-8o E 

278.305 

102.1 

U 
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The emission gap of 15 km s 1 width at ?’lsr +15 kms 1 seen in the SMA HCN J— 3-2 data 
in Figure[2]is likely to be an artifact resulting from the limited uv coverage in our observations. The 
single-dish 13 CO map shows that the emission in this velocity range extends beyond the primary 
beam of the SMA, implying that it is resolved out in the synthesis observation. 


3.2. Double-lobed Structure: Tight Anti-correlation between the Broad Velocity 

Lobes 


The double-lobed structure detected with the ASTE 10-m telescope observation is also clearly 
recognizable in the SMA image with the improved resolution. The P-V diagram of HCN J— 3-2 
presented in Figure [2] shows that the two broad velocity lobes are sharply separated in space at 
l = —0°.306, with a narrow overlapping ^lsr range from +15 kms -1 to +35 kms -1 . Figure [3] shows 
the peak flux density maps of the HCN J— 3-2 for the 'T’lsr ranges of the blue-shifted lobe (from 
—35 kms -1 to +15 kms -1 ), the intermediate velocity range (from +15 kms -1 to +40 kms -1 ), 
and the red-shifted lobe (from +40 kms -1 to +110 kms -1 ), along with the single-dish HCN J— 4-3 
maps for the same velocity ranges. The SMA images reveal that the two broad velocity lobes 
have a ridge-like shape; each broad velocity lobe has a vertical structure near the field center, as 
schematically shown in the upper rightmost panel. Hereafter we refer to this pair of parallel ridges 
as the central ridges. The emissions from the red-shifted lobe primarily arise from the central ridge 
part. The blue-shifted lobe also shows bright emissions in the eastern region behind the central 
ridge, making the entire lobe appear inverse-T’ shaped. The central clump is located toward the 
mid-point of the blue-shifted central ridge. The emissions from the intermediate velocity range 
(—15 kms -1 to +40 kms -1 ) are mostly extensions from the two broad velocity lobes, and no 
notable spatial structures intrinsic to this velocity range are detected. 


Figure 0] shows composite color images constructed from the two broad velocity lobes, where 
the 3 a emission regions of the blue- and red-shifted lobes are shown in cyan and red, respectively. 
As reported in Tanaka et al. & the two velocity lobes in the single-dish image show a good 
spatial anti-correlation, with a narrow overlapping region of an about half the beam size. Now we 
recognize that the anti-correlation is much tighter in the interferometric images; the central ridges 
correspond exactly to the overlapping region in the single-dish image, but they have almost no 
spatial gaps and only few overlaps between each other. Thus the two broad velocity lobes appear 
to contact each other at a well-defined interface on the plane of the sky. A natural explanation for 
this remarkable structure is that the two components with different velocities physically interact 
with each other, that is, they collide. Meanwhile, our results do not completely exclude other 
interpretations, as the paired ridges also seem to compose one elongated ring. We will discuss this 
issue further in later sections mm 14X21) . 


The emissions of the central ridges partly mix in the intermediate velocity range. Figure [5] 
shows the P-V diagrams of the SMA HCN J— 3-2 data along horizontal paths across the northern 
end (b = — 0°.59), midpoint (b — — 0°.64), and southern end (b = — 0°.69) of the ridges. We find 
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several features bridging the velocity gap between the central ridges along these paths. Along the 
central strip, the two velocity components are smoothly connected by emission with a velocity 
gradient in a 'T’lsr range from —30 kms -1 to +60 kms -1 , except for the artificial emission gap 
noted in the previous subsection. A bridging feature is also detected along the northern strip in a 
narrower v^sr range from +20 kms -1 to +40 kms -1 . The southern strip lacks detectable emission 
in the intermediate velocity range, although it is possibly because of the missing flux. 


3.3. Broad HCN Emission 

Figure [ 6 ] shows the HCN J— 3-2 P-Fdiagrams along the Galactic latitude at every 10" Galactic 
longitude, with contours of CH 3 OH 9_i-8o E emission overlaid, where the HCN and CH 3 OH spectra 
are smoothed with Gaussian kernels of 2 kms -1 and 5 kms -1 FWHM, respectively. The broad 
velocity lobes are resolved into a number of smaller broad emissions. The red-shifted lobe (+lsr > 
+40 kms -1 ) is found to comprise at least 6 broad features with about 10" sizes and velocity widths 
of 40 - 50 kms -1 in FWZI, the systemic velocities of which are distributed in a 'Ulsr range from 
+35 kms -1 to +90 kms -1 without apparent order. The blue-shifted lobe (+lsr +15 kms -1 ) does 
not seem to have prominent broad features but rather seems to be an ensemble of relatively narrow 
clumps, although this is at least partly because of the artificial emission gap at vlsr ~ +15 kms -1 ; 
any broad features in this velocity range would be unrecognizable as they would be divided by this 
artificial gap. The narrow emissions in the 'L’lsr range from —30 kms -1 to —20 kms -1 could be 
the low-velocity ends of such obscured broad features. 

The broad emission features are predominantly associated with the central ridges mentioned 
in the previous subsection, which correspond to the Galactic longitude range from —0°.303 to 
—0°.308 in Figure [ 6 l They are rather randomly distributed in the position-velocity space, in which 
we cannot detect definitive signs of systematic motion over the entire clump, such as coherent 
expansion, rotation, and linear collimated flow. 

We label 8 bright methanol peaks as A-H in Figure [ 6 j the positions of which on the l-b plane 
are also shown in Figure 01 Peak A is the methanol peak in the central clump. All of the methanol 
peaks have velocity widths of 10-20 kms -1 , which are narrower than those of the HCN emissions. 
They are primarily associated with the broad HCN features in the central ridges, except for peak 
B, which is located well outside the ridges, and peak H, which is located inside the red-shifted ridge 
but has an inconsistent blue-shifted velocity. We list the positions of the methanol peaks in Table El 
along with their center LSR velocities and FWHM velocity widths determined by Gaussian fitting 
of the 52 ~ 4 i E spectra. 


3.4. Central Clump 


Figure [7^ shows the SMA images of the central clump in the continuum, and the N 2 H + J— 3-2 
and CH 3 OH 9o~8i E flux densities integrated over a velocity range from —10 kms -1 to +20 kms -1 . 
Table [3] summarizes the properties of the central clump. The clump has approximately the same 
elliptic shape in the continuum and in the N 2 H + maps; they are elongated in the northeast- 
southwest direction with FWHM sizes of 13" x 8 ". The point source-like methanol 9o~8i E emission 
is offset from the center position of the clump toward the south-western edge, i.e., toward the 
boundary between the two velocity components. 


We estimate the dust luminosity mass of the central clump by referring to the Herschel archival 
data. First, we determine the dust temperature (Ta) by analyzing the spectral energy distribution 
(SED) constructed from the 70-500 /im flux densities taken from the Herschel data archive. Figure 
E b shows the SED for the central 40" square region, where contributions from extended emission 
unassociated with the cloud are subtracted by taking the background position at ±60" offset from 
the cloud center in the Galactic longitude and latitude. The error bars include the calibration 
errors in the Herschel data and uncertainties of the background levels estimated using the standard 
deviations among the flux densities toward the four reference positions. The SED is consistently 
fitted by a single-component gray-body model with a dust temperature (Tj) of 18 K and the 
dust emissivity index (/?) assumed to be 2.0. Although a simultaneous SED fitting of T& and (3 
suggests a higher f3 value of 2.8 , we fix 0 at the conservative value of 2.0. Th is result is consistent 
with former observations (e.g. 


Pierce-Price et al. 


2000 ; 


Molinari et al, 


2011 ) indicating that the 


dust temperatures are uniformly ^ 20 K throughout the CMZ. No significant excess in short 
wavelengths is detected. The mass of the central clump is estimated to be 5.8 x 10 2 M o from this 
dust temperature and the 1.1 mm continuum luminosity (0.38 Jy) measured for the entire clump 
with our SMA data. We obtain the mean density no = 5 x 10 5 cm -3 from the mass and size of the 
clump by roughly assuming the volume to be ^ 7 rr 3 . This is cons istent with the density estimate 
obtained from the excitation analysis of the HCN intensities ( Tanaka et al. 20141 ). 


Table 2 . Methanol Peaks 



1 

(°) 

b 

H 

^LSR 
(km s -1 ) 

Av 

(km s -1 ) 

Reference 

A 

-0.305 

-0.064 

1.1 ±0.2 

12.8 ±0.3 

Central Clump, 36 GHz maser candidate 1 

B 

-0.300 

-0.058 

0.1 ±0.4 

12.1 ±0.9 


C 

-0.309 

-0.071 

45.4 ± 1.0 

22.0 ±2.4 


D 

-0.308 

-0.065 

26.0 ± 1.0 

23.7 ±2.5 


E 

-0.309 

-0.068 

54.0 ± 1.7 

17.7 ±4.5 


F 

-0.309 

-0.062 

59.5 ±0.5 

6.8 ± 1.2 


G 

-0.309 

-0.060 

45.4 ±0.7 

13.6 ± 1.8 


H 

-0.309 

-0.065 

-0.2 ±0.9 

17.8 ±2.3 



^Yusef-Zadeh et al.l (l2013h 
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We estimate the viral mass (My t) of t he centr al clump by using the formula Myx/M© = 


210(r/pc) • (dp/kms -1 ) 2 used in Miyazaki fe Tsubo i (2000) for the CS clumps in the CMZ. The 
radius r is estimated as \D^/6 in a j • 0 m i n = 0.15 pc, where D — 8.3 kpc is the distance to the source 
and 0 ma j, min = 10".2, 6".6 are the major and minor axis lengths corrected for the instrumental 
broadening due to the finite beam size, respectively. The FWHM velocity width dv is measured to 
be 12.4 kms -1 from the averaged N 2 H + spectrum, after correction for broadening due to unresolved 
hyperfine structures. These values yield Myx = 4.8 x 10 3 M©. The virial parameter <a v i r = Myx/M 
is 8.3, which is only slightly smaller than the typical value for the CMZ clouds, ~ 10, estimated 
from the CS J=l-0 observations by Miyazaki fc Tsuboi (2000b By considering the rather large 
uncertainties in dust and gas properties, we conclude that the central clump is not significantly 
different from the typical CMZ clouds in <a v i r . 


3.5. Methanol Peaks 

Table 0] lists the fluxes of the lines detected toward the methanol peaks A-H. Toward the 
methanol peaks outside the central clump, i.e., peaks B-H, only two methanol transitions, 52-4 1 E 
and 9_i-8oE, and the N 2 H + line are detected above the 3 -a levels. We use a stacking analysis 
to estimate the fluxes of the undetected lines toward peaks B-H; their spectra are co-added after 

Table 3. Properties of the Central Clump and the Methanol Peaks 


Parameter 

Value 

Entire Clump 

Dust temperature (T^) 

18.4 ±0.4 K 

Radius (r) 1 ’ 2 

0.15 ±0.01 pc 

FWHM velocity width (An) 1,3 

12.4 ±0.2 km s -1 

Dust luminosity mass (M) 

5.8 x 10 2 M© 

Virial mass (My t) 

4.8 x 10 3 M© 

Mean density 

5 X 10 5 cm -3 

Mehtanol Peak-A 

Methanol rotation temperature 

28 ± 1 K 

Gas density 4 

1Q6.5-7.5 cm -3 

Mehtanol Peaks-B-H 

Methanol rotation temperature 

28 ± 14 K 


^rom gaussian fitting of the N2H+ image. 

2 Corrected for the beam broadening. 

3 Corrected for the broadening by unresolved hyperfine 
structures. 

4 By assuming the gas kinetic temperature to be 40-100 


K. 
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correction for the difference in the central velocities. All the lines except 13 CS J— 5-4 are detected 
above the 3 -a levels in the stacked spectra, the fluxes of which are listed in the bottom line of Table 

m 


The left column of Figure [8] shows methanol rotation diagrams constructed from the spectra 
toward the methanol peak A and the spectra stacked for peaks B-H, in which the upper state 
column densities per unit velocity (d7V u /di;) divided by the statistic weights (g u ) are plotted against 
their excitation energies for the four methanol transitions. The dN u /dv values are calculated 
from averaged flux densities over the three central velocity channels around the 52-4 1 E peaks by 
assuming the optically thin condition. For the central clump, the flux densities of three out of the 
four transitions are successfully fitted by a single-temperature local thermodynamical equilibrium 
(LTE) model with a rotation temperature (T rot ) of 28 K, but the remaining 9 _i-8q E transition has 


an inconsistently high value possibly because the 9 _i-8qE transition is masing. lYanagida et al. 


E3 detected class-I maser emissions of this methanol line that are associated with molecular 
outflow sources in the infrared dark cloud G34.43+00.24 MM3. Peak A also coin cides with a source 
in the list of 36 GHz class-I methanol maser candidates by Yusef-Zadeh et al. (2 0131 ). Therefore, 
we exclude the 9_i-8o E transition in our data from the excitation analysis, although it is yet to be 
confirmed whether the transition is actually masing because its spectral profile is similar to those of 
thermal lines and lacks sharp peaks typical for maser profiles. For the stacked data of the methanol 
peaks B-H, a single temperature LTE fit is less successful, and the uncertainty in the rotational 
temperature is accordingly larger. The best fit T rot value is again 28 K. 


We perform a non-LTE excitation analysis to estimate the physical conditions of the methanol 
peaks by assuming the optically thin condition. The rate coefficients for the radi ative and collisional 
transitions are taken from the Leiden atomic and molecular database (jSchoier et al .1 12 OO 5 I ). The 
right column of Figure [El shows the ranges of gas kinetic temperature (T^m) and molecular hydrogen 
density (nn 2 ) that reproduce ±1 a values around the measured flux density ratios among the 6j- 
52 E, 52~4i E, and 9 q-8i E transitions. In this analysis, we conservatively assume 10% relative errors 


Table 4. Line Fluxes toward Methanol Peaks in Jy kms 1 



61-62 E 

CH3OH 

52—4iE 9 q— 81 E 

9-i-8oE 

n 2 h+ 

3-2 

o-H 2 CS 

8l,7-7l,6 

OCS 

23-22 

13 cs 

6-5 

A 

10.6 ±1.4 

60.4 ± 1.3 

9.1 ± 1.3 

51.0 ± 1.3 

125.6 ± 1.5 

14.6 ±1.0 

9.8 ± 1.6 

9.6 ±1.4 

B 


17.5 ±1.1 


12.2 ± 1.3 

11.4 ± 1.6 




C 


16.7 ± 1.6 


15.9 ± 1.9 

6.2 ±0.9 




D 


22.7 ±2.0 


19.7 ± 1.9 

11.3 ± 1.3 




E 


7.8 ±1.6 


10.0 ± 1.9 

9.4 ±2.0 




F 


6.2 ±0.9 


10.0 ± 1.1 

8.7 ± 1.3 




G 


10.1 ± 1.1 


8.8 ± 1.2 

9.5 ± 1.3 




H 


15.8 ±1.7 


14.9 ± 1.3 

4.1 ± 1.2 




B-H (stacked) 

1.3 ±0.3 

13.8 ±0.6 

2.0 ±0.3 

12.8 ± 0.6 

8.5 ±0.8 

3.9 ±0.5 

2.3 ±0.4 
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in the flux densities (i.e. 14% in the ratios) in addition to the spectral noises. In the result for peak 
A, the 6 i-52 E/52-4 i E and 9o~8i E/52-4i E ratios give a consistent parameter range, showing that 

h a reasonable assumption for the gas temperature, 40 < T < 100 K 


- 10 6 - 5-7 - 5 cm ' 3 


(Martin et al. 

2004; 

Nagai et al. 

2007) 


Nagai et al. 2007). For the stacked data for the methanol peaks B-H, the 


parameter ranges calculated from the two flux density ratios do not overlap. If we adopt the result 
for the 9o~8i E/52-4 1 E ratio, which has a higher S/N ratio than the other ratio, we obtain a nn 2 
range similar to that deduced for peak A. The parameters derived in this subsection are listed in 
Table [3j 


DISCUSSION 


4.1. Kinematics of CO—0.30 — 0.07 


In our previous paper (Tanaka et al. 2014|), we presented four hypotheses for the origin of the 
broad velocity emission of CO—0.30—0.07: cloud-cloud collision, expansion driven by blast wave, 
molecular outflow, and rotation around an invisible massive object. We discuss these hypotheses 
in this subsection. 


4-1.1. Cloud-Cloud Collision 


With the interferometric HCN J— 3-2 image taken with the SMA, we found remarkably tight 
anti-correlation between the two broad velocity lobes comprising CO—0.30—0.07, which is much 
tighter than that previously found with single-dish ASTE observations; the characteristic broad 
velocity emissions are predominantly confined in the central ridges running through the cloud 
center, and the individual ridges in the two velocity components are tightly anti-correlated with 
each other on the plane of the sky. Thus the entire cloud appears to be sharply divided into the 
two velocity components by the linear boundary defined by the central ridges. 


We can give a natural explanation for this structure by assuming that_the cloud is a pair of 
colliding clumps. In the numerical model of colliding clouds by Gilden (119841 ). planar compressed 
layers appear behind the shock fronts propagating away from the collision interface and into the 
individual clouds. The central ridges of CO—0.30—0.07, which run along the interface between the 
two velocity components, are consistent with this model. The tu rbulent velocity in the post-shock 
region is expected to be enhanced up to the shock velocity (Inoue et al. 2013|). For CO—0.30—0.07, 
the minimum value of the collision velocity v co \ is roughly estimated to be 20 kms -1 , i.e., half of 
the 'T’lsr difference between the two velocity components, although the edge-on geometry of the 
collision interface may indicate that the tangential collision velocity is faster than the line-of-sight 
velocity. The velocity dispersion of 18 kms -1 inside the red-shifted lobe is consistent with this 
collision velocity within the uncertainty in the tangential velocity. 
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We can identify a spatially extended part of each of the colliding pair in the single-dish 13 CO 
J— 3-2 map. Figure [9b is the P-V diagram of the 13 C0 J— 3-2 data along the diagonal path across 
the cloud indicated in Figure [9b- The figure shows multiple components in the velocity range 
from —50 kms -1 to +90 kms -1 in addition to the broad emissions of CO—0.30—0.07. The blue- 
shifted lobe of CO—0.30—0.07 is obviously associated with the + 15-kms _1 cloud, the most intense 
and most widespread component in the map. The red-shifted lobe also accompanies a spatially 
extended, narrow component at vlsr ~ +55 kms -1 . This +55-kms -1 cloud is elongated from the 
south-western direction and is truncated at the position where it intersects the red-shifted lobe. 
Figure [9b shows a composite color image constructed from the single-dish 13 CO data, in which the 
+15- and +55-kms -1 clouds are shown in blue and red, respectively. The figure illustrates that the 
broad velocity lobes observed in the SMA data (indicated by the contours in the figure) are well 
correlated with the boundary between the +15- and +55-kms -1 clouds. This configuration, i.e., 
compact broad emissions at the contact_points_of two_clouds with narrower line widths, is typical 
for a cloud-cloud collision system (e.g. Duarte-Cabral et al. 2011). 


A major argument against the cloud-cloud collision scenario is that a part of the shocked 
region is actually extended well behind the interface between the two clouds in our data; the blue- 
shifted lobe also has bright emissions in the eastern part of the cloud, where it lacks a red-shifted 
counterpart (Figure [3j ). In order to explain this feature, we have to adopt additional, ad hoc 
assumptions such as that the red-shifted counterpart is somehow invisible or has already dissipated 
or that there is another source of turbulence in the eastern side of the +15-kms -1 cloud. 


4.1.2. Expanding Molecular Ring Model 

An expanding ring driven by a blast wave is also a good model for the cloud kinematics, if 
we ignore the spatial anti-correlation between the central ridges by regarding it as coincidental. In 
the ring model, the two ridges are assumed to be portions of an ellipse elongated in the Galactic 
north-south direction, as indicated in Figure [T0h. The figure also shows the first velocity moment 
map of the SMA HCN J— 3-2 data. The cloud has an overall velocity gradient across the major 
axis of the ring, which is interpreted as expansion rather than rotation by adopting the natural 
assumption that the ring is nearly circular in the three-dimensional space with a large inclination 
from the plane of the sky. The other possible kinematics, that is, a model of a highly elliptical 
rotating ring inclined around the minor axis, is discussed in the next sub-subsection (I4.1.3|) . 

Figure [TOb shows the angular variation of the first moment of the ring defined in Figure flOh. 
along with the model curve for an expanding ring with an inclination of 62°.5 and expansion velocity 
of 45 kms -1 . Although we have failed to detect clear signs of coherent expansion motion in the 
P-V diagrams (Figures [2] and [6j), the first moment variation in Figure \10\ is successfully fitted 
by the expansion model; both the model and observation show steep gradients at the tangential 
points (0^0°, 180°) and rather flat curves in the other parts. The failure in the detection of 
coherent expansion in the three-dimensional P-V images could be explained as being the result 
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of disturbance caused by interaction with non-uniformly distributed outer media; such a situation 
would be expected for an expanding ring driven by explosion in the vicinity of a dense molecular 
cloud. 


An advantage of the expanding-ring model over the collision model is that it can easily explain 
the feature unassociated with the central ridges without additional assumptions because an isotropic 
explosion can create shocked regions also outside the expanding ring. A problem with the expansion 
model is that the region inside the assumed ring is actually not an emission cavity; it shows low- 
level emissions smoothly bridging the velocity gap between the central ridges, as shown in Figure 
[5l The central ridges are most smoothly bridged at their midpoints by a horizontal structure, 
making the ring appear to be shaped like the symbol 0. The expanding-ring model cannot provide 
a simple explanation for the horizontal bar structure. In addition, the absence of visible energy 
sources capable of driving the fast expansion motion, the kineticenergy of which is estimated to be 
4 x 10 49 ergs from the total mass of CO—0.30—0.07 (10 3 M 0 ; Tanaka et ah 20141 ) and the expansion 
velocity (45 kms -1 ), remains as a major problem; t he cloud is not associa ted with a ny sources in 


the published catalogs for radio (INord et al. 2004; Lazio & Cordes 2008) or X-ray ( Muno et al, 


2006) sources. We also note that the first moment variation displayed in Figure flOb is insufficient 


to distinguish between the collision model and expanding-ring model because the observed pattern 
is expected for any cases involving two distinct velocity components of similar sizes along the path. 


4-1.3. Other Scenarios 


The other two models in llTanaka et al.ll (2014), namely the rotating-ring/disk and the hyper- 
energetic outflow scenarios, do not fit the observation. If we assume that the central ridges represent 
an elliptic Keplarian orbit around a compact massive object, the large eccentricity implies that 
the attractor is very close to either end of the major axis and that very steep velocity gradients 
should appear at the pericenter passage as a result. In Figure flOh we show two rotation models 
corresponding to the two possible attractor loci (i.e., the northern and southern focal points) 
with various inclinations and a fixed maximum line-of-sight velocity of 80 kms -1 . Neither model 
reproduces the observed line of sight velocities, especially those near the pericenter passage. The 
models with low inclinations, which show milder velocity variation near the pericenter, may seem 
to fit the observation better, but they require unrealistically large attractor masses: 1 x 10 6 Mq 
and 1 x 10 7 M 0 for inclinations of 30° and 10°, respectively. 

The bipolar-outflow model is also unsuccessful; the broad emissions are not linearly aligned as 
expected from the outflow model. 
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4-1.4- The Origin of the Broad Emissions 

We have shown that the cloud-cloud collision model gives the most reasonable explanation 
for the cloud kinematics and the absence of visible energy sources with fewer assumption than the 
other three models. We can safely rule out the rotating-ring and bipolar-outflow models with our 
observations. The expanding-ring model cannot be completely excluded, because it reproduces the 
observed velocity structure of the cloud very well for a limited spatial range. Further investigation 
of the shock structure with a better resolution will help to determine whether the shocked gas 
represents a collision interface or is shaped as an elongated expanding ring. 


4.2. Implications of the Shock-Triggered Star Formation 


Our SMA observations toward CO—0.30—0.07 have revealed a remarkable structure quite 
different from that known from previous observations: a compact, quiescent clump with bright 
N 2 H+ and continuum emissions co-exists with the highly turbulent gas traced by the HCN and 
HCO + lines. We also detected at least 8 point source-like, bright methanol peaks, including one in 
the central clump. The densities of the methanol peaks are deduced to be 10 6 ' 5 7,5 cm 3 from the 
excitation analysis with the multi-transition methanol data. The high density, compactness in size 
(< 0.1 pc in radius), and bright methanol emissions suggest their similarity to typical star-forming 
cores. Hence, we consider that CO—0.30—0.07 may provide an ideal test case for star formation in 
the CMZ, w here stars generally form under highly turbule nt environments relative to the Galactic 
disk region (jKruijssen et al.ll2014l ; iRathborne et al.ll2014bl ). 


The central clump and all of the methanol peaks are associated with broad HCN emissions, 
especially those aligned along the boundary between the two velocity components (Figured]). This 
suggests that they are density enhancements created via compression by the shock formed at the 
boundary, irrespective of whether the origin of the shock is cloud-cloud collision or blast-wave 
impact. This assumption is corroborated by the widths of the methanol lines; they are approxi- 


mately 2-3 times that exp 

ected from their sizes 

5 and the size-linewidth relationship for the CMZ 

(Mivazaki & Tsuboi 

2000 

Oka et al. 2001). The theory of turbulence-regulated star formation 

(Mac Low & Klessen 

2004; 

Krumholz & McKee 

2005; 

Padoan &; Nordlund 

2011 

) suggests that the 


density enhancements in turbule nt c louds start collapsing once they become self-gravitating and 
eventually form stars. Inoue & Fukui (2013) show that high turbulent pressure and enhanced mag¬ 
netic field strength in cloud-cloud collision systems set a preferable environment for the formation 
of massive stars. We may expect such a triggered star formation process in CO—0.30—0.07. 


As indicated by the high a v [ T of ^ 10, the entire central clump is not self-gravitating but likely 
confined by external ram pressure. It is unclear whether the central clump and methanol peaks 
have any self-gravitating substructures capable of forming stars. Here we investigate whether the 
methanol peak A in the central clump is sufficiently dense to become self-gravitating, according to 






































15 


the modeling in Federrath & Klessen (2012). The critical density n cr ;t is given by 

n C rit/«o ~ A ■ a v ; r • A4 2 ■ (l + /3 -1 ) -1 , 


(1) 


where no and Ai are the mean density of the cloud and r.m.s. sonic Mach number of the turbulence, 
respectively, and A is a constant close to unity. The value f3 = 2c^/c\ is the thermal-to-magnetic 
pressure ratio, where c s and c\ are sound speed and Alfven velocity, respectively. If the magnetic 
field is negligible (/3 —>► oo), the critical density is 8 x 10 8 cm -3 for the central clump when - m = 
50 K. We obtain a lower critical density of ^ 2 x 10 7 cm -3 by assuming that ca = 3 kms -1 from 
the typical values of the density (nn 2 = 10 3 ' 5 cm -3 ; INagai et al. 2007) and magnetic field strength 
(B 0.1 mG: iFerrierel 120091 . and references therein) in the CMZ clouds and that the magnetic field 


is frozen-in to the gas. The latter value is comparable to the density deduced for peak A, implying 
that peak A is possibly dense enough to be self-gravitating, although we will need spatially resolved 
observations toward the methanol peaks in order to make a more conclusive evaluation. Here, we 
note that the strong tidal field in the CMZ, which is not taken into account in Equation (PQ), does 
not significantly change the above discussion. The critical density for self-gravity to o vercome the 
tidal shear is approximately 10 4 cm -3 at a Galactocentric radius of 100 pc (jYusef-Zadeh et al 


2007), indicating that the tidal shear is negligible in the density range considered. 


It may seem that stars are already forming in the methanol peaks in CO—0.30—0.07 because 
of the richness of the molecular species detected toward them; bright methanol emission is a major 
in dicato r of m assive star formation, and OCS and H 2 CS are also classified as hot-core-type molecules 
(llmmer et al. 2014). However, we refrain from identifying the methanol peaks as hot cores because 
the hot-core-like chemi stry is not limited t o star -forming regions in the CMZ; the hot core chemistry 
without hot cores (IRequena-Torres et al . 2006) is rather widespread throughout the region, where 
the molecules formed on the grain surfaces are efficiently evaporated into the gas p ha se by co smic- 


rays and/or shocks in addition to heating by young massive stars fjReqnena- Torres et al.l 12006 


Menten et al . 2008; Yusef-Zadeh et al. 2013)- We estimate the fractional methanol abundance to 
molecular hydrogen toward peak A to be 2 x 10 -7 on the basis of the methanol and continuum flux 
densities by assuming an ortho/para ratio of unity and the LTE condition with a T rot of 28 K. This 
high abundance is typical for hot cores but is also well within the range of the abundances measured 
toward non-star-forming clouds in the CMZ, suggesting that the high methanol abundance toward 
the methanol peaks might simply reflect the typical chemical composition in the CMZ rather than 
indicate hot-core formation. The class-I methanol maser candidates in the 9_i-8o E transition 
could be further evidence for massive star formation, but it is also possible that they are excited 
directly by the shock impact rather than by outflow from youn g ste llar objects. In fact, the 36 GHz 
class-I methanol maser candidates identified in Yusef-Zadeh et al. (2013) show poor correlation 
with known star-forming regions. Internal heating sources are also undetected in the dust SED of 
the central core (Figure Eb). 

The collision-induced cluster form ation scenario is a popular explanation for the formation of 


super star clusters (SSCs) in the CMZ ( Hasegawa et al.lll994l : iHiguchi et al.ll2014l ). CO—0.30—0.07 


is a new candidate for the site of triggered star formation at its early stage, although it is evident 
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that the cloud, the mass of which is ^ 10 3 M©, is only capable of forming a cluster of at most a few 
100 M©. We note that the situation of CO—0.30—0.07 is similar to that of GO.253+0.016, or the 
“brick”, the infrared dark clo ud in the CMZ often re f erred to as a progeni t or of a SSC similar to th e 


Arches and the Qui ntuplet (jLongmore et al.1 l2012al : iHiguchi et al.l 12014 : iRathborne et al.1 l2014al ) . 
Higuchi et al. (2014) show sulfur monoxide (SO) cores with large velocity dispersion (a v ) aligned 


along an arc-shaped structure, which they argue is an imprint of past impact by a smaller cloud. 
The methanol peaks confined in a thin shocked layer in CO—0.30—0.07 could be smaller versions 
of these large-cr^ SO cores. 


5. SUMMARY 

By utilizing the SMA, we have performed synthesis mapping toward the HVCC CO—0.30—0.07, 
a peculiar molecular cloud with an extremely broad velocity width in the CMZ. We summarize the 
main results below. 

1. The cloud has two distinctive components: a highly turbulent, spatially extended emission 
traced by HCN and HCO + and a central clump of 0.15 pc radius with bright N 2 EU and 
continuum emissions. The central clump has a narrow velocity width of 12 kms -1 in FWHM, 
which is in striking contrast to the extremely broad emission of the extended component 
spanning over a 145 kms -1 velocity range. 

2. The extended component has a pair of parallel vertical ridges running through the center of 
the cloud. These central ridges are tightly anti-correlated with each other in both space and 
velocity, thereby dividing the entire cloud into two components with line-of-sight velocities 
differing by 40 kms -1 . 

3. The marked broad emissions of the cloud primarily arise from the central ridges, which are 
again resolved into a number of smaller broad emission features with widths of 40-50 kms -1 
in the HCN image. These broad HCN emission features are rather randomly distributed in 
the position-velocity space without apparent order. No pattern indicating systematic motion 
such as coherent expansion, rotation, or bipolar outflow is detected. 

4. The mass of the central clump is estimated to be 5.8 x 10 2 M© from the continuum luminosity. 
The virial parameter a v [ T is not significantly different from those of typical CMZ clouds, ^10. 

5. At least 8 bright methanol peaks are detected toward the central clump and in the cen¬ 
tral ridges, accompanying several other hot-core-type molecular lines. Their methanol rota¬ 
tion temperatures are estimated to be 28 K under the LTE assumption. Non-LTE analysis 
for the methanol peak of the central clump shows that the molecular hydrogen density is 
1q6.5— 7.5 cm -3 unc [ er the assumption that T^[ n = 40-100 K. 
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6. CO—0.30—0.07 is located at a position where two larger clouds at ^lsr — +15 and +55 
kms -1 contact each other, and the central ridges with broad emissions are well correlated 
with the boundary between them. The tight anti-correlation between the ridges suggest that 
CO—0.30—0.07 could be a site of cloud-cloud collision. Meanwhile, another model in which 
the central ridges are assumed to be parts of a molecular ring expanding at a velocity of 45 
kms -1 cannot be excluded, as it explains the angular variation of the line-of-sight velocity 
well. The other two models proposed in lTanaka et al. ( 2014 ). namely the bipolar-outflow and 
the rotating-ring models, do not fit the observations. 


7. The central clump and the methanol peaks, which have small radii (< 0.1 pc), densities 
marginally high enough to be self-gravitating, and rich abundances of hot-core-type molecules, 
are similar to typical star-forming cores in the Galactic disk region. These objects are predom¬ 
inantly aligned along the central ridges and hence are likely to be features created by shock 
impact. We speculate that the cloud is at an early phase of shock-triggered star formation, 
although conclusive evidence for star formation activities is yet to be obtained. 


The authors are grateful to an anonymous reviewer for his/her helpful comments. We also 
thank the SMA staff for their support during the observations reported here. The Submillimeter 
Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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Fig. 1.— Maps of the peak flux densities of the observed lines and the continuum flux densities. The 
peak flux densities are calculated from spectra binned over 10 velocity channels. The synthesized 
beam is represented by a filled ellipse on the HCN J— 3-2 panel. 
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HCN J=3-2 N 2 H + J = 3-2 
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Fig. 2.— Galactic longitude-velocity diagram of CO—0.30—0.07 averaged over the full latitude 
range in HCN J— 3-2 and N 2 H + J— 3-2 obtained with the SMA observations, along with the HCN 
J=4-3 and 13 CO J— 3-2 images obtained with the ASTE observations. The dark hatched region on 
the leftmost panel indicates the velocity range affected by the artificial velocity gap in the synthesis 
observation. 
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Fig. 3.— Peak flux density/brightness temperature maps for the blue-shifted lobe (leftmost; 'T’lsr 
range from —35 kms -1 to +15kms _1 ), intermediate velocity range (second from left; from +15 
kms -1 to +40 kms -1 ), and red-shifted lobe (third from left; from +40 kms -1 to +110 kms -1 ). 
The top and bottom rows are for the SMA HCN J— 3-2 and ASTE J=4-3 data, respectively. The 
rightmost panel of the upper row is a schematic diagram showing the positions of the central ridges 
in the blue- and red-shifted lobes (shown in blue and red, respectively) and the central clump. The 
coverage of the SMA observation is superposed on the HCN J=4-3 maps. 
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Fig. 4.— Composite color images of the two broad velocity lobes, in which the 3 a emission regions 
of the blue-shifted and red-shifted lobes are shown in cyan and red, respectively, and overlapping 
regions are shown in white. The region corresponding to the central ridges is encircled by yellow 
dashed lines. The methanol peaks A-H are indicated by cross marks. 
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Galactic Longitude (deg) 

Fig. 5.— HCN J— 3-2 position-velocity diagrams along the horizontal strips passing thorough 
northern end (b = — 0°.059), midpoint (—0°.064), and southern end (—0°.069) of the central ridges. 
The bridging emission between the two broad velocity lobes are indicated by dashed lines. The 
artificial velocity gap in the synthesis observation is hatched in gray. 
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red-shifted lobe 



Fig. 6 .— Galactic latitude-^LSR diagrams of HCN J— 3-2 at every 10" Galactic longitude, with 
contours of CH 3 OH Jk = 9_i-8o E overlaid. The contour levels are (3, 5,10, and 20) x a (a = 0.10 
Jybeam -1 ). The HCN and CH 3 OH spectra are smoothed with Gaussian kernels of 2 kms -1 and 5 
kms -1 FWHM, respectively. Bright methanol peaks are labeled as A-H. The hatched regions are 
for the artificial velocity gap in the synthesis observation. 
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Fig. 7.— (a) SMA images of the central clump in the continuum flux density (left) and the 
N 2 H+ J— 3-2 velocity-integrated flux density (right). The integration range is from —10 kms -1 to 
+20 kms -1 . The overlaid contours are the CH 3 OH Jk = 9o~8i E flux integrated over the same 
velocity range, drawn at every 1 a from 5 a (a = 1.5 Jybeam -1 kms -1 ). The boundary between 
the central ridges is schematically drawn with a dashed line on each panel, (b) Spectral energy 
distribution of the central clump in the far-infrared wavelengths, constructed from 70 - 500 fi m flux 
density from the Herschel data archive. The best-fit gray-body curve (solid line) and 1 -a deviation 
(hatched region) are overlaid. 
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a) Methanol peak A (Central Clump) 



b) Mathanol peaks B-H 



Fig. 8.— (a) Methanol rotation diagram toward peak A, with the best-fit curve of a single¬ 
temperature LTE analysis overlaid (left), and the results of non-LTE excitation analysis, in which 
the parameter ranges calculated from the observed flux density ratios of Jk — 6 i-52 E and 9o~8i E 
to 52~4 i E are hatched in red and gray, respectively (right). In the latter analysis, we assume 10% 
relative error for each flux density in addition to the spectral noise levels, (b) Same as (a) but for 
the stacked data for the methanol peaks B-H. 
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Fig. 9.— (a) Position-^LSR diagram of the ASTE 13 CO J— 3-2 data along the path defined in 
panel b, with contours drawn at every 2 a starting from 3 a (a = 0.11 K). The velocity ranges 
of the +15-kms -1 and +55-kms -1 clouds are indicated, (b) Composite color image constructed 
from the 13 CO J— 3-2 integrated intensities of the +55-kms -1 and +15-kms -1 clouds (red and 
blue color channels, respectively). Contours of the peak HCN J— 3-2 flux densities of the red- 
and blue-shifted lobes in the SMA data are overlaid and shown in red and blue, respectively. The 
contours are drawn at 3 cr intervals starting from 7 a (a = 0.28 Jybeam -1 ). 
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Fig. 10.— a) The first moment map of the SMA HCN J— 3-2 data, with the colors and brightness 
representing the first velocity moment and peak flux density, respectively, (b) Angular variation of 
the first moment along the ellipse defined in panel a. The definition of the position angle 6 is also 
given in panel a. The overlaid curve is for an expanding-ring model with an expansion velocity v e ^ p 
of 45 kms -1 and inclination angle of 62°.5. (c) Same as (b) but with model curves for rotating-ring 
models with various assumptions for the inclination angle i. The maximum line-of-sight velocity is 
fixed at 80 kms -1 . The attractor loci are assumed to be the northern and southern focal points in 
models-1 (left) and -2 (right), respectively. 













